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Abstract. The ACDM cosmological model assumes the existence of a small cosmological 
constant in order to explain the observed accelerating cosmic expansion. Despite the dramatic 
improvement of the quality of cosmological data during the last decade it remains the simplest 
model that fits remarkably well (almost) all cosmological observations. In this talk I review the 
increasingly successful fits provided by ACDM on recent geometric probe data of the cosmic 
expansion. I also briefly discuss some emerging shortcomings of the model in attempting to 
fit specific classes of data (eg cosmic velocity dipole flows and cluster halo profiles). Finally, I 
summarize recent results on the theoretically predicted matter overdensity (Sm = ^^) evolution 
(a dynamical probe of the cosmic expansion), emphasizing its scale and gauge dependence 
on large cosmological scales in the context of general relativity. A new scale dependent 
parametrization which describes accurately the growth rate of perturbations even on scales 
larger than lQQh~^ Mpc is shown to be a straightforward generalization of the well known scale 
independent parametrization /(a) = S7m(ci)^ valid on smaller cosmological scales. 



1. Introduction 

Converging geometrical [H El El HI El E] and dynamical O El El UHl E] cosmological observations 
indicate that the universe has entered a phase of accelerating expansion. This expansion could 
be driven either by the negative pressure of a homogeneous dark energy component |12j or by a 
modified gravitational force which is repulsive on large cosmological scales. Alternatively it could 
simply be due to an unusually large underdensity (void) on scales of about IGpc which induces 
an apparent isotropic accelerating expansion to observers located approximately at the center 
of the void [13]. All of the above three classes of models attempting to explain the observed 
accelerating expansion suffer from fine tuning problems which are reflected in the unexpected 
values of the parameters required to fit the cosmological expansion data. For example dark 
energy originating from an evolving scalar field requires a scalar field mass of order lQ~^^GeV 
which is several orders of magnitude smaller from expectations based on particle physics theories. 
Similarly, models assuming the existence of large voids require fine tuning of the location of the 
earth based observer within about 2QMpc from the IGpc large void while the formation of such 
a large void is also very unlikely in the context of standard cosmology |13j . 

The simplest and most successful model consistent with (almost) all cosmological data is the 
ACDM model |14J. This model is based on a modified version of Einstein's equations of the form 



(1) 



where A is a new free parameter, the cosmological constant. When the new term Ag^;^ is placed 
on the left hand side (Ihs) of equation ([1]) it is interpreted as a modification of the gravitational 
law corresponding to an effective Newton potential [13 



V{r) = -^-Ar' (2) 

where — Ar^ corresponds to a new repulsive gravitational term. When the new term Ag^^ is 
placed on the right hand side (rhs) of equation ([1]) it is interpreted as a new contribution to 
the energy momentum tensor T^i, corresponding to dark energy with constant density pA = 
and constant negative pressure p\ = —p\. 

Consistency with cosmological data requires a fine tuned value for p\ 

p^^^'^ ~ {lf)-^'^GeVf ~ 2 X IQ-^^erg/cm^ (3) 

A physically motivated origin of the cosmological constant could be the energy of the quantum 
field vacuum which is predicted to have a constant diverging energy density and negative 
pressure. Under the assumption of a proper cutoff the quantum vacuum can be made finite 
and play the role of a cosmological constant. A natural scale for this cutoff is the Planck scale 
leading to a vacuum density 

p^P^ ~ {IQ^^GeVf ~ 2 X 10^i°er5/cm3 (4) 

which is 120 orders of magnitude larger than the observed value. 

Despite of the unnaturally small value of the observed p^^^^ when compared to the anticipated 

value /ojf'^ in the context of vacuum energy, the cosmological constant of equation ([TJ has some 
unique physically motivated features. In particular the Ihs of equation ([1]) is the most general 
second rank tensor which is 

• local 

• coordinate covariant 

• divergenceless (needed for energy momentum conservation) 

• symmetric 

An additional important attractive feature of ACDM is simplicity. It is the only model 
based on General Relativity (GR) which (assuming flatness) involves a single free parameter 
(JIa = ^ ~ ^ ~ ^Om where pc is the present day critical density required for a flat universe) 
and its predicted expansion rate H{z) = ^(z) as a function of redshift z 

H{zf = Hi [Oom(l + zf + (1 - Oom)] (5) 

is currently consistent with all cosmological expansion probes (geometric and dynamical). It 
is therefore clear that the identification of potential confiicts of ACDM with cosmological data 
is a prerequisite before seriously considering alternative more complicated models unless such 
models are free from any type of fine tuning. Unfortunately no such model is currently known. 
Cosmological observations testing the ACDM model include two classes of probes 

• Geometric probes which measure directly the cosmic metric (eg Type la supernovae as 
standard candles or Cosmic Microwave Background (CMB) spectrum peaks and Baryon 
Acoustic Oscillations (BAO) which use the last scattering horizon scale as a standard ruler) 



• Dynamical probes which measure simultaneously the cosmic metric and the gravitational 
law on cosmological scales. The main dynamical probe is the evolution of linear dark matter 
cosmological perturbations 5m{k,z) = ^^{k,z) as a function of scale and redshift. This 
evolution can be probed either directly through weak lensing observations [15] or indirectly 
though the large scale power spectrum of luminus matter at various redshifts |16j . 

Geometric probes provide currently the most accurate information about the cosmic expansion 
rate H{z) as a function of redshift. In the first part of this brief review (section 2) I discuss the 
following three basic questions 

(i) What is the figure of merit |17j (constraining power) of various recent Type la supernovae 
(Snia) datasets and how does it compare with the corresponding figure of standard ruler 
(CMB+BAO) geometric probe data? 

(ii) What is the consistency level of various geometric probe datasets with ACDM? 

(iii) What is the level of consistency between standard candle SnIa datasets and standard ruler 
CMB+BAO geometric probes? 

A few cosmological observations which appear to be in some tension [T8] with the predictions of 
ACDM will also be discussed in section 2. The most interesting of these observations appears 
to be the observed dipole velocity flows on scales larger than 50h^^ Mpc which appear to be a 
factor of about four larger than the ACDM model predictions assuming normalization of the 
matter fluctuation power spectrum using the WMAP5 CMB data[19l [20l [21]. Other puzzling 
observations include the high redshift brightness of SnIa \22l , the cluster halo profiles [231 121] 
and the emptiness of voids [251 IMl 127] . 

The theoretically predicted linear evolution of cosmological perturbations 5m{k,z) as a 
function of redshift and scale, depends sensitively on both the expansion rate and the 
gravitational law on cosmological scales. This dual sensitivity makes the linear evolution of 
density perturbations a particularly useful dynamical cosmological probe. Even though current 
observational constraints on 5m{k,z) are not as powerful as geometric constraints on H{z) this 
is expected to change dramatically in the next decade [28j. 

The time evolution of matter density perturbations 6m{k,t) is described by the well known 
equation (see eg [8]) 

5„ + 2H6„^ - iTTGpmfik, t)5„, = (6) 

where in the context of general relativity f{k,t) = 1 and 6m{k,t) = 5m{t) becomes independent 
of the scale k. In fact, this scale independence of 6m is occasionally considered to be a signature 
of validity of general relativity [29]. However, the derivation of equation (jlSp is based on 
two important assumptions whose validity is questionable on scales (~ 200h~^ Mpc). These 
assumptions are the following: 

(i) The scale of the perturbation 6m is significantly smaller than the Hubble scale at all times 
during the perturbation evolution. 

(ii) The choice of gauge plays a minor role in the form of the evolution equation for 6m' 

An important question discussed in the second part of this review (section 3) is the following: 
What is the level of validity of these assumptions on large cosmological scales? It can be 
shown that the validity of both assumptions breaks down rapidly as the scale increases beyond 
200h-^Mpc. 

2. ACDM confronts Recent Geometric Cosmological Data 

A useful way to test the ACDM model is to consider a generalized parametrization of the 
dark energy equation of state which includes ACDM as a special parameter case and find the 



likelihood of the parameters corresponding to ACDM in this context. A commonly used such 
parametrization is the Chevallier-Polarski, Linder (CPL)[30l 131] ansatz 

2; 

w = wq + wi{l - a) = Wo + wi- (7) 

± ~\~ z 

which reduces to ACDM for {wq,wi) = (—1,0). In this section, I show the ranking of the 
six latest Type la supernova (Snia) datasets (see Table [T]) (Constitution (C), Union (U), 
ESSENCE (Davis) (E), Gold06 (G), SNLS lyr (S) and SDSS-II (D)) in the context of the CPL 
parametrization ([7]) according to their Figure of Merit (FoM)[17j. their consistency with the 
cosmological constant (ACDM) and their consistency with standard rulers (Cosmic Microwave 
Background (CMB) and Baryon Acoustic Oscillations (BAO)). The datasets considered are 
shown in Table ([T|) along with some useful features such as the redshift range or the subsets of 
each set [32]. 

Assuming a CPL parametrization for 'w{z) (equation ([7])) it is possible to apply the maximum 
likelihood method separately for standard rulers (CMB+BAO) and standard candles (SnIa) 
assuming flatness. The corresponding late time form of H{z) for the CPL parametrization is 

H\z) = Hi[no^{l + zf + 

+ (l-f]om)(l + ^)'^'+"°+"^)eWir]. (8) 



Table 1. The datasets used in the present analysis. See respective references for details on the 
sources of the SnIa data points. 



Dataset 



Date Released Redshift Range # of SnIa Filtered subsets included 



SNLSl L33j 

Goidoe m\ 



ESSENCE [36], [37] 



Union 
Constitution [3] 

SDSS 



2005 
2006 

2007 

2008 
2009 

2009 



0.015 < z < 1.01 115 

0.024 <z< 1.76 182 

0.016 <z< 1.76 192 

0.015 <z< 1.55 307 

0.015 <z< 1.55 397 

0.022 < z< 1.55 288 



SNLS [33], LR 
SNLSl [3J], 

HST [35J, SCP d] 
HZSST [2] 



SNLSl [33J, HST 
ESSENCEt36j,t37] 



GoldOe [35], ESSENCEf36], 

m 

Union [38], CfA3[3] 

Nearby |l0], SDSS-II [39j, 
ESSENCE [36], SNLS [33], 
HST 135] 



The SnIa observations use light curve fitters [40] to provide the apparent magnitude m{z) of 
the supernovae at peak brightness. The resulting apparent magnitude m{z) is related to the 
dimensionless luminosity distance Dl{z) through [HI [32] 



mth{z) = M{M,Ho) + 5logio{DL{z)) , 



(9) 



where 

DUz) = {l + z) r dz'-—^ ^ (10) 

is the Hubble free luminosity distance (Hodi), and M is the magnitude zero point offset and 
depends on the absolute magnitude M and on the present Hubble parameter Hq as 

M = M + 5logio{-^) + 2b = 
Mpc 

= M - 5/oc/io/i + 42.38. (11) 

The parameter M is the absolute magnitude which is assumed to be constant after proper 
corrections (using light curve fitters) have been implemented in m{z). 

The theoretical model parameters are determined by using the maximum likelihood method 
ie by minimizing the quantity 

2 ir, \ if^obs{zi) - Hthizi))^ , . 

XSnlai^, Wo, Wi) = 2^ 2 (12) 

i=l i 

where is the number of Snia of the dataset and cj^ j are the errors due to flux uncertainties, 
intrinsic dispersion of SnIa absolute magnitude and peculiar velocity dispersion. These errors 
are assumed to be Gaussian and uncorrelated. The theoretical distance modulus is defined as 

IJ-thizi) = mth{zi) - M = Uogio{DL{z)) + jiQ , (13) 

where 

/io = 42.38 - 5/0510 /i, (14) 

The steps followed for the usual minimization of ()12p in terms of its parameters are described 
in detail in [32l El ISl |l2] . 

In the context of constraints from standard rulers from CMB spectrum peaks |41l [52] , the 
maximum hkelihood method is used and the datapoints (/?, /a, lOOOfo/i^) of Ref. g] (WMAP5) 
where i?, la are two shift parameters [H]. In the case of BAO, the maximum likelihood method 
is also applied [H] using the datapoints of Ref. [13] (SDSS5). For comparison, the more recent 
data of Ref. [S] (SDSS7) have also been considered with only minor differences in the results 
(slightly reduced consistency with ACDM in the context of the CPL parametrization but no 
change in the ranking sequences). 

The Figure of Merit (FoM) is a useful measure of the effectiveness of a set of data 
in constraining cosmological parameters. In the case of two parameters (as for the CPL 
parametrization) it is defined as the reciprocal area of the 95.4% contour (see Fig. [2]), in 
parameter space {wq^wi) [17J. Clearly, the larger the FoM the more effective the dataset in 
constraining the parameters {wq^wi). Fig. [1] shows the FoM in terms of the number of the SnIa 
data for the datasets of Table [TJ Clearly, the FoM is an increasing function of the number of 
SnIa in the datasets. An exception to this rule is the ESSENCE dataset which has a slightly 
smaller FoM compared to the GoldOG dataset even though it has a larger number of SnIa. A 
possible origin of this effect is the fact that the FoM does not depend only on the total number 
of SnIa of the dataset but mainly on the number of SnIa at low and high redshifts (the redshift 
space distribution of the ESSENCE data includes more data at intermediate redshifts than the 
GoldOG dataset). 

In order to study the consistency of the various SnIa datasets with the cosmological constant 
and the standard rulers we consider the distance in units of a (cj-distance d^) of the best fit 
point to a model with parameters {wq^wi)^ where this reference point can be either ACDM 
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Figure 1. The Figure of Merit (FoM) in terms of the number of the SnIa data for riom = 0.28 
using the CPL parametrization. 



or some other reference point, (see Fig. The u-distance can be found by converting 

Ax^ = - xLn to i-e- solving 01] 

1 - r(l, AxV2)/r(l) = Erf (djV2) (15) 

for da (cr-distance), where Ax^ is the difference between the best-fit and the reference point 
{wo,wi) (eg ACDM) and Erf() is the error function. The right hand side of Eq. ()15p comes 

from integrating f ""J!" — 7= e 2^ dx, where n is the desired number of as, while the left hand 

side corresponds to the Cumulative Distribution Function (CDF) of a distribution [H] with 
two degrees of freedom. Note that Eq. p5|) is only valid for the two parameters {wi^^w\) and 
should be generalized accordingly for more parameters [Uj. In the special case of n = 1 or 
n = 2 we obtain the well known results Axi^- = 2.30 and Ax2o- = 6.18 valid for two parameter 
parametrizations |44j . 

Even though the cr-distance is not a commonly used statistic it is quite useful because it 
can directly give information about probability of a given region in parameter space. The 
integer values of sigma distance (Icr and 2cr) are commonly used to draw the corresponding 
contours in parameter space. The extension of this statistic to non-integer values is used to find 
the specific contours that go through particular reference points of parameter space and thus 
estimate quantitatively the consistency of these points. The advantage of using the cr-distance 
instead of Ax^ is the fact that the a-distance takes into account the number of parameters of 
the parameterizations and can therefore be directly translated into probability for each point in 
parameter space. This is not possible for Ax^ because it does not include information about 
the number of parameters of the parameterizations considered. 

It is straightforward to apply the likelihood method to find the trajectory of the best 
fit point {wq^w\) in parameter space as Oom varies in the range Oom £ [0.2,0.36]. These 
trajectories obtained for each of the datasets of Table [1] and also for the standard ruler CMB- 
BAO (WMAP5+SDSS5, WMAP5-hSDSS7) data are shown in Fig. El These trajectories can 
not be used to directly rank the datasets according to their consistency with any given reference 
point in parameter space (e.g. ACDM) because they contain no information about the 68.3% 
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b igure 2. The 68.3%(lcj) - 95.4% (2cj) confidence contours in the wq — wi plane based on 
parametrization ^ for the ESSENCE (left) and SNLSl datasets (right) for Oom = 0.24. The 
arrows indicate the cj-distance of ACDM (green points: {wq, wi) = (—1, 0)) to the best fit points 
(red points). 
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Figure 3. Trajectories of the best fit points {wo,wi) obtained for each of the datasets of Table 
Hand also for the standard ruler CMB-BAO (WMAP5-FSDSS5 and WMAP5-FSDSS7) data as 
r^Om varies in the range r^om G [0.2,0.36]. The arrows in the best fit lines indicate the direction 
of growing r^om- Note that for the SDSS5 data the standard ruler best fit parameters stretch 
out to {wq,wi) ~ (2, —30) for ilom — 0.36, whereas for the SDSS7 data {wo,wi) ~ (0.90, —20). 



(Icr) contours. However, they provide useful hints for the trend of the best fit parameters as 
riom varies. For example, such a trend is the increase of the best fit value of the slope wi as the 
prior of decreases towards the value 0.2 or that the best fit value of wq remains less than 
— 1 for all datasets except of the Constitution compilation. 



■"o.20 0.25 0.30 0.35 °i,.2o 0.25 0.30 0.35 

inO il ..,() 

Figure 4. a: cr-distances (io-i(riom; — 1, 0) off the ACDM point (reference point) from the best fit of 
each dataset in Table [T] Notice that the best fit of each dataset, except for the Gold06 and the SDSS-II 
(MLCS2k2), are minimized at similar values of fiom- The cr-distance d^{flo^; —1, 0) between the standard 
ruler best fits and ACDM are also shown as a function of 51om (black and grey solid lines), b: Similar 
to (a) for the dynamical dark energy reference point (uiqjWi) — (—1.4,2). Notice that the a-distances 
daii^Om] —1-4, 2) are minimized at more widely separated values of riom- 

The ranking sequence of the datasets of Table [T] with respect to any reference point in 
parameter space can be studied quantitatively using the u-distance statistic discussed above. In 
order to test the sensitivity of the ranking sequence of datasets with respect to the choice of 
consistency reference point {wojWi) we consider two such reference points: {wq^wi) = (—1,0) 
(ACDM) and (wq, wi) = (—1.4, 2) which corresponds to dynamical dark energy with a w{z) that 
crosses the line w = —1. It should be noted that there is nothing special about the parameter 
point (-1.4,2). We have selected it as a representative of a wide region in parameter space (upper 
left from ACDM) which corresponds to dynamical dark energy crossing the phantom divide line 
w = —1. Any other point in the same parameter region would lead to similar results and the 
same ranking of datasets. This particular parameter region is interesting because it is spanned 
by the best fit trajectories and it also mildly favored by the Gold06 dataset (see Fig. E]). 

The resulting da{^om) for each dataset of Table [1] are shown in Figs. H] in the range 
Oom G [0.2,0.36]. 

Clearly, there are values of fiom that minimize the cr-distance dai^oia) between the best 
fit of each dataset and the reference point. These values of Oom maximize the consistency 
of the datasets with the given reference point in this range of r^om- The minima cr-distances 
dai^Om) for each dataset, corresponding to maximum consistency with ACDM along with the 
corresponding values of Oom are shown (properly ranked) in Table [21 The corresponding results 
for the reference point {wo,wi) = (—1.4,2) are shown in Table [3l 

The following comments can be made with respect to the results shown in Figs. [Hand in the 
corresponding Tables [2 [3l 

(i) The consistency with ACDM of all datasets, except GoldOG and SDSS-II when using the 
MLCS2K2 method, is maximized in a narrow range of Qom S [0.26, 0.29] which also includes 
the value of fiom favored by standard rulers. On the other hand, the consistency with the 
dynamical dark energy point (wqjWi) = (—1.4,2) is maximized over a wider range of Oom 
(Slom £ [0.27,0.35]) thus decreasing the consistency among the datasets in the context of 
dynamical dark energy. 

(ii) The ranking sequence changes dramatically when the consistency with the dynamical dark 
energy is considered as a reference point instead of ACDM (Table[3|). Essentially the ranking 
is reversed! Thus, the choice of the consistency reference point plays an important role in 
determining the ranking sequence of the datasets (see also Fig. [J]). 



Table 2. Minimum cr-distance (i™"(r2oJ^"; —1,0) from the best fit point for each of the datasets 
to the ACDM point. Also listed are the corresponding values of fiom) and the best fit parameters 
(wq, Wl) (see also Fig. [3^). The SDSS-II (MLCS2K2) data showed no minimum of da with respect 
to riom in the range ilom G [0.2,0.36]. We thus have simply displayed the lowest value of in 
the corresponding range of r^om- 



Dataset 






Wq 


Wl 


SNLSl 


0.004 


0.260 


-1.03 


0.16 


SDSS-II (SALT2) 


0.084 


0.270 


-1.09 


0.51 


Constitution 


0.114 


0.285 


-0.91 


-0.54 


ESSENCE 


0.227 


0.270 


-1.20 


1.04 


Union 


0.525 


0.285 


-1.25 


1.40 


SDSS-II (MLCS2K2) 


0.623 


0.360 


-1.06 


0.93 


GoldOe 


0.950 


0.345 


-1.56 


2.80 


CMB+BAO (SDSS5) 
CMB+BAO (SDSS7) 


0.200 
0.588 


0.272 
0.272 


-1.15 
-1.30 


0.51 
0.97 



Table 3. The minimum cr-distances d™" to the reference point {wq^wi) = (—1.4,2), the values 
of riom at which the minimum distance is attained, and the best fit parameters {wq,wi) at 
r^o^" are displayed for each dataset (see also Fig. Hja). We omit the rows corresponding to 
CMB-I-BAO data as the resulting u-distance is always ^ 1 due to the dominance of the dark 
energy at early times. The SDSS-II (MLCS2K2) data showed no minimum of d^ with respect 
to r^om in the range ilom G [0.2,0.36]. We thus have simply displayed the lowest value of d^ in 
the corresponding range of $7om- 



Dataset 


"cr 




Wq 


Wl 


GoldOe 


0.11 


0.345 


-1.56 


2.80 


Union 


0.17 


0.300 


-1.26 


1.25 


ESSENCE 


0.19 


0.275 


-1.21 


0.99 


SNLSl 


0.42 


0.290 


-1.04 


-0.26 


SDSS-II (SALT2) 


0.63 


0.275 


-1.10 


0.46 


SDSS-II (MLCS2K2) 


0.87 


0.360 


-1.06 


0.93 


Constitution 


1.77 


0.315 


-0.88 


-1.32 



(iii) The SDSS-II dataset obtained with the MLCS2k2 fitter has some peculiar features compared 
to other datasets. In particular it favors particularly high values of riom (f^Om — 0.4) while 
for r^om < 0.3 its consistency with ACDM is significantly reduced to a level of 3(T or larger 
{da > 3). In addition, the trajectory of its best fit parameter point as rJom varies is 
perpendicular to the corresponding trajectory of most other datasets (see Fig. H]). 

It is straightforward to apply the cr-distance statistic to rank the Snia datasets according to 
their consistency with standard ruler CMB-BAO data. We simply use as a consistency reference 
point the best fit point {wq, wi)^^ for standard rulers obtained as described in section 2 using the 
WMAP5-I-SDSS5 data. In this case, the location of the reference point {wq, wi)^^ in parameter 




Figure 5. a-distance do-i{^om', {wo,wi) ) between the reference point {wo,wi) of standard 
rulers (using the WMAP5+SDSS5 data) and the best fit of each dataset as a function of 
for the datasets of Table [TJ Using the SDSS7 data, the minimum distances are found in the 
range Qoni e [0.29,0.31]. 



space depends on fiom but this does not complicate the analysis. The (T-distance between the 
reference point {wo,wi)^^ and the best fit of each dataset is shown in Fig. [5] as a function of 
for the datasets of Table [TJ These distances are minimized for values of riom that are different 
for each dataset but they are all in the narrow range r^om £ [0.27,0.3]. 

These minimum distances along with the corresponding value of Qom are shown in Table [J] for 
each dataset (properly ranked according to consistency with standard rulers). Notice that the 
ranking sequence for the consistency with standard rulers is practically identical to the ranking 
sequence of the consistency with ACDM (Table [2]) but differs from the ranking sequence of the 
consistency with dynamical dark energy (Table [3]). This is an interesting feature of the data in 
favor of ACDM. 



Table 4. Consistency with standard rulers. Minimum cj-distance (i™"(f]om; (^O) "U^i)"^^) 
between best fit parameters for each dataset, {wqjWi), and best fit for standard rulers, 
{wo,wi)^^. ^^"'(ilom; {'Wo,wi)^^) is minimized at r^g^" (see Fig. E]). 



Dataset 




ntmin 


Wo 


Wl 




^SR 


SNLSl 


0.003 


0.280 


-1.04 


-0.10 


-1.07 


0.08 


SDSS-II (SALT2) 


0.058 


0.280 


-1.11 


0.40 


-1.07 


0.08 


ESSENCE 


0.087 


0.275 


-1.21 


0.99 


-1.12 


0.38 


Constitution 


0.121 


0.295 


-0.90 


-0.76 


-0.84 


-1.28 


Union 


0.681 


0.280 


-1.24 


1.44 


-1.07 


0.08 


GoldOe 


1.976 


0.280 


-1.38 


2.75 


-1.07 


0.08 


SDSS-II{MLCS2K2) 


3.342 


0.285 


-0.92 


1.02 


-1.00 


-0.28 



It is therefore clear that the ACDM cosmological model is a well defined, simple and predictive 
model which is consistent with the majority of current cosmological observations. Despite of 
these successes there are specific cosmological observations which differ from the predictions of 
ACDM at a level of 2a or higher. The observations conflicting the WMAP5 normalized ACDM 
model at a level of 2a or larger include the following [18j : 

• Large Scale Velocity Flows (ACDM predicts significantly smaller amplitude and scale of 
flows than what observations indicate) [19\ I21j . The probability of consistency with 
ACDM is about 1%. 

• Brightness of Type la Supernovae (Snia) at High Redshift z (ACDM predicts fainter Snia 
at High 2;) [22]. The probability of consistency with ACDM is about 3 — 5% for the Union 
and Gold06 datasets. 

• Emptiness of Voids (ACDM predicts more dwarf or irregular galaxies in voids than 
observed) [251 ESI EZ]. 

• Profiles of Cluster Haloes (ACDM predicts shallow low concentration and density profiles 
in contrast to observations which indicate denser high concentration cluster haloes) |23[l24j. 

• Profiles of Galaxy Haloes [35] (ACDM predicts halo mass profiles with cuspy cores and low 
outer density while lensing and dynamical observations indicate a central core of constant 
density and a flatfish high dark mass density outer profile), 

• Sizable Population of Disk Galaxies [46] (ACDM predicts a smaller fraction of disk galaxies 
due to recent mergers expected to disrupt cold rotationally supported disks). 

Even though some of the puzzles discussed here may be resolved by more complete observations 
or astrophysical effects, the possible requirement of more fundamental modifications of the 
ACDM model remains valid. 



3. Dynamical Probes: The Grow^th Function 

3.1. Growth of Perturbations in General Relativity: Beyond the sub-Hubble Approximation 
As discussed in the introduction, the evolution of the matter overdensity 6m{k,a) (the growth 
function) consists a useful probe of both the expansion rate and the gravitational law on large 
scales. The standard parametrization of the linear growth function 6m{k, a) is scale independent 
and is obtained by introducing a growth index 7 defined through the growth rate /(a) by 

/o(a) = = (16) 

a in a 

where a = j^:^ is the scale factor and 

^^mW- ^^^^2 [in 

is the ratio of the matter density to the critical density when the universe has scale-factor a where 
Hq is Hubble constant and ilom is ratio of mass density to critical density. This parametrization 
[17] provides an excellent fit to the evolution equation for ^^(a) in general relativity and in the 
small scale (sub-Hubble) approximation 

5m + 2H6m - ^TTGpmSm = (18) 

where an overdot denotes the derivative with respect to time and pm is the matter density. 
Changing variables from t to In a we obtain the evolution equation for the growth factor / as 

f' + f + fi§2+^) = l^rn (19) 



where ' = d/dlna. For dark energy models in a flat universe with a slowly varying equation of 
state w{a) = = wq, the solution of eq. (fT9]) is well approximated by eq. (fT6]) with [17] 

(20) 

which reduces to 7 = ^ for the ACDM case {wq = —1). It is therefore clear that the 
observational determination of the growth index 7 can be used to test ACDM [9]. It has 
been shown [11] that even in the context of dynamical dark energy models consistent with Type 
la supernovae (Snia) observations the parameter 7 does not vary by more than 5% from its 
ACDM value. However, in the context of modified gravity models 7 can vary by as much as 
30% (e.g. for the DGP model[48j ^dgp — 0.68 [11]) while scale dependence is also usually 
introduced[lHlEDl[ID]- 

Current observational constraints on 7 are based on redshift distortions of galaxy power 
spectra [JT], the rms mass fluctuation a%{z) inferred from galaxy and Ly — a surveys at 
various redshifts [SSJ [S3] , weak lensing statistics [31] , baryon acoustic oscillations [5S] , X-ray 
luminous galaxy clusters [SB], Integrated Sachs- Wolfe (ISW) effect [ST] etc. Unfortunately, the 
currently available data are limited in number and accuracy and come mainly from the first two 
categories. They involve significant error bars and non-trivial assumptions that hinder a reliable 
determination of 7. Thus, the current constraints on 7 are fairly weak [9] and are expressed as 

7 = 0.674^0:195 (21) 

This however is expected to change in the next few years when more detailed weak lensing 
surveys [15] are anticipated to narrow significantly the above range [28]. 

A crucial assumption made in the derivation of eq. ()18p in the context of general relativistic 
metric perturbations is the assumption that the scale of the perturbations is significantly smaller 
than the Hubble scale [37]. This assumption however does not lead to a good approximation 
on scales larger than about 100/i~^-/Vfpc [58] I59j . In order to demonstrate this fact consider the 
perturbed metric of spacetime which takes the form (in the Newtonian gauge): 

ds'^ = -(1 + 2^)dt^ + (1 - 2^)a^-fijdx'dx^ , (22) 

where 7jj is the metric of the spatial section and we are ignoring anisotropic stresses. The 
evolution of density perturbations on all scales is dictated by combining the background 
equations 

= -^{Pm + Pde) (23) 

p = -3H{p + p) (24) 

(assuming a flat universe with only pressureless dark matter and (non-clustering) dark energy) 
with the perturbed linear order Einstein equations in the Newtonian gauge [6OJ {pm and pde are 
the matter and dark energy densities respectively while p = wp is the pressure). The resulting 
(anisotropic stress-free) equations are of the form 

$ = -4H^ + 8TrGpdeWde'^ (25) 

6 = 3<^> + ^Vf (26) 
Vf = (27) 



with constraint equations 



h2 

m{H^ + ^) + ^^ = -ATrGSpm (28) 
{H<^ + ^) = -AnGp^n^Vf (29) 

where $ is the Newtonian potential, Vf = —va {v is the velocity potential for dark matter). 
Clearly, equations (p5]) - (j27|) involve a scale k dependence in contrast to the small scale 
approximate equation (jlSp which is scale independent. 

The derivation of equation (jlSp in the context of general relativity is made using the sub- 
Hubble approximation. The linear matter overdensity 6pm may be expressed [60] in terms of 
the gravitational potential $ and the background variables as follows: 

- 4TTG5pm = + 3^/2$ + 3H<^> (30) 



a2 



In the sub- Hubble (small scale) approximation (^ ^ H"^) equation (pOj) takes the form 

- AitG5p^ = (31) 

where a slowly varying gravitational potential <I> has also been assumed. 

The general relativistic equations (|25p - ()29p lead to the following equation for the matter 
overdensity 5m 

6m + 2HSm + = (32) 

which also expresses the conservation of the perturbed energy momentum tensor for matter. 
Using the sub-Hubble approximation ()3ip we obtain the scale independent approximate equation 
(jlSp . On the other hand, if we avoid this approximation in equation ()30p . solve for <1> (ignoring 
the time derivative) and substitute in equation ([32]) we obtain the following scale dependent 
evolution equation for 5m [SSI EH] : 

5m + 2H5m - = (33) 



where 



1 + C{a,k) 



?,a^H{af 



aa,k) = ^ (34) 



The solution of equation ()33p provides a much better approximation to the full linear general 
relativistic system ()25p - ()27p up to horizon scales. On scales larger than the horizon even equation 
(|33p breaks down since on these scales, the time derivative of $ can not be ignored. 

Given the successful approximation of the solution of ()33p to the exact linear general 
relativistic solution, it becomes important to construct a scale dependent parametrization that 
is analogous to ()16p and solves (approximately) ()33p for all scales k. In order to construct such 
a parametrization we focus on the matter dominated era when most of the growth occurs and 
express ^(a, fc) as 

aa,k) = '-^ (35) 
Equation §3^ may be expressed in terms of the growth factor / — "^dina^ the form 

/' + /+(2-|0..w)/ = |^f|gL. (30) 
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Figure 6. a: The growth rate / obtained from the solution of the general relativistic 
system (k = O.Olh^^ Mpc, r^om = 0.3, ACDM , dotted line) compared with the 
scale independent parametrization (continuous line) and the corresponding generalized scale 
dependent parametrization (thick dashed line) b: Similar to a. for k = 0.004/i~^Mpc. c: 
Similar to a. for k = O.OOlh^^ Mpc. 

where ' = and we have assumed ACDM for H{a). 

For sub-Hubble scales ■^(/c, a) — )• and equation ([36]) reduces to (fT9|) whose solution is well 
approximated by ()16p with 7 = It is straightforward to show using an expansion method 
|59] that the parametrization 



is an approximate solution of equation (j36p and provides a good approximation to the solution 
of the general relativistic system (j25p - (j27p up to horizon scales. 

The accuracy of the parameterization p7p is demonstrated in Fig. [6] where we compare the 
form of the scale dependent parametrization Eq. ()37p with the general relativistic numerical 
solution and with the standard parametrization for three different scales. It is clear that up to 
approximately the Hubble scale {k ~ 0.001hMpc~^) the scale dependent parametrization ()37p 
is accurate at a level better than 5% at least up a to redshift z = 10. It can be shown that this 
parametrization is also a good approximation to the growth rate not only in the ACDM model 
but also in the case of dynamically evolving dark energy. Thus the parametrization described 
by equation ()37p approximates very well the full linear general relativistic solution up to horizon 
scales and differs significantly from the standard parametrization on scales larger than about 
WOh-^Mpc. Its use may play an important role when comparing cosmological structure data 
on large scales with theoretical model predictions. 

As shown in Fig. [6] the range of scales where the standard scale independent parametrization 
()16p starts breaking down involves scales larger than WOh~^ Mpc {k < 0.01hMpc~^). This scale 
is much smaller than the present Hubble scale which is about 3000/i~ ^Mpc. Therefore, the sub- 
Hubble approximation breaks down at much smaller scales than the anticipated Hubble scale. 
This is due to the fact that at early times during the matter era when most growth occurs, 
the comoving Hubble scale is significantly smaller compared to its present value. In fact, at 
recombination it is close to 200h~^ Mpc. Therefore, scales that remain at the sub-Hubble level 
during the whole time when the growth of matter perturbations take place are only scales below 
200h-^Mpc. 

The generalized Poisson equation (|30p in the Newtonian gauge may be solved for the slowly 
varying $ leading to a generalized gravitational potential in Fourier space of the form 

^Poisson f ],\ 

i + P 



where ^P°^^^°^ ^j-^'j [g the usual smah scale gravitational potential which emerges as a solution of 
the usual Poisson equation. In coordinate space the generalized large scale potential potential 
has a Yukawa form with a Hubble scale cutoff namely 

$G«(r) = _^e-v^^^- (39) 
or 

It is this general relativistic potential that should be compared with corresponding modified 
gravity potentials on scales comparable to the Hubble scale. 

In a modified gravity theory the Poisson equation on large scales is modified for two reasons: 
First due to the modified gravitational law and second due to the Hubble scale effects that are 
present also in the case of general relativity as discussed above. The result is a Poisson equation 
of the form 

= -ATTGpmSrnik, a)f{k, o) (40) 

where the scale dependent function f{k, a) incorporates the scale dependence due to both Hubble 
scale effects and scale dependence due to modification of the gravitational law. This function is 
to be compared with the scale dependence due to pure Hubble scale effects fcRik, a) = fc) 
present in general relativity. Therefore a signature of modified gravity would not be a scale 
dependence of the evolution of matter density perturbations as claimed occasionally in the 
liter ature|29j. but a deviation from the scale dependence predicted on large scale in the context 
of general relativity. 



3.2. Gauge Dependence 

The scale dependence of the matter density growth rate obtained above was based on calculations 
made in the Newtonian gauge. This is a physically important gauge because it corresponds to a 
time slicing of isotropic expansion. Nevertheless since the matter density perturbation 6mik,t) 
is a gauge dependent quantity, it is important to clarify to what extend do the results of section 
3.1 persist in a gauge different from the Newtonian gauge. 

Another important gauge is the synchronous gauge which corresponds to a time slicing 
obtained by the matter local rest frame everywhere in space (the free falling observer frame). 
The synchronous gauge is generally considered to be the most efficient reference system for 
doing calculations. It is used in many modern cosmology codes for calculating the evolution of 
cosmological perturbations (eg CMBFAST [61]). The line element of the perturbed spacetime 
in the synchronous gauge is given by 

ds^ = a\T){-dT^ + ((5y- + h,j)dx'dx^) (41) 

where r is the conformal time. It is straightforward to derive the growth equation for 
<Jm = 5pm/ Pm in a matter dominated universe in the synchronous gauge to obtain [601 EH] 

¥^ + 2H5^^ - A^p^GS^^ = (42) 



This growth equation is exact in the synchronous gauge in the case of matter domination and 
involves no scale dependence as in the case of equation ([33]) of the Newtonian gauge. This scale 
independence is an artifact of the particular time slicing of the synchronous gauge which is a 
good approximation on small scales but is unable to capture the horizon scale effects modifying 
the growth function on large scales. 

Nevertheless equations ()33p and ()42p clearly agree on small scales where ^ — )• 0. Therefore, for 
larger scales {k < O.OlhMpc) the question that arises is the following: What is the proper gauge 
to use when comparing with observations? This question was recently addressed in Ref. |62] 



where a gauge invariant observable replacement was obtained for 6m- This observable 6obsik,t) 
involves the matter density perturbation 6m{k, t) corrected for redshift distortions due to peculiar 
velocities and gravitational potential. It also includes volume and position corrections. The final 
expression however is complicated and makes the theoretical predictions based on it not easy to 
implement and manipulate. 

Alternatively, the gauge- invariant (GI) approach to the cosmological perturbations evolution, 
pioneered by Bardeen|63J may be used to identify observables (e.g., Refs.fMl [55]). The most 
general form of the line element for a spatially flat background and scalar metric perturbations 
can be written as |65] 

ds'^ = a^{{l + 2$)dr2 - 2B\idx'dT - [5ij - 2{^5ij - E\^^)]dx'dx^ }, (43) 

where a and r are the conformal cosmic expansion scale factor and the conformal cosmic time; "|" 
denotes the background three-dimensional covariant derivative. The corresponding perturbed 
energy-momentum tensor Tjf has the form 

Tq" = PmU\i , 

T\ = -p^(U-B)\i, 

= -pm^^J , (44) 

where pm is the unperturbed pressureless matter density; U and S determine velocity 
perturbation and anisotropic shear perturbation. 

A gauge-invariant matter density perturbation may be constructed as [631 El] 

6g^s ^6m + 3-{U-B) , (45) 
a 

S^^^ coincides with the density perturbation dm"^^^ in the Comoving Time-orthogonal Gauge 
(CTG, in which U = B = 0), which denotes the density perturbation relative to the spacelike 
hypersurface which represents the matter local rest frame everywhere [63]. This quantity also 
coincides with the density perturbation 6^ in the Synchronous Gauge (in which "if = ^ = B = 
0) for the pressureless matter system. In other words, S^^^ denotes the density perturbation 
relative to the observers everywhere comoving with the matter. These free falling observers 
do not experience the isotropic expanding background of the universe because the peculiar 
velocity of matter is distinct from the Hubble flow. Thus S'^^ has physical significance only for 
perturbations on scales small compared to the Hubble scale. 

An alternative gauge-invariant variable is more closely related to the matter overdensity in 
the Newtonian gauge. This gauge-invariant perturbation variable is of the form |63 [ 1641 [65] . 

6g'^ ^6m + -{B-E) = 6m- S-{B - E) . (46) 
p a 

and has important advantages over . coincides with the density perturbation in 

the Newtonian Gauge (NG), in which B = E = 0. 

In addition to the gauge invariant perturbation variable it is straightforward to construct two 
gauge-invariant scalar potentials and ijj, both of which become the same as the gravitational 
potential in the Newtonian limit. These are constructed from metric perturbations as follows 

m- 

<P^'^--{B-E) , 
a 

^ = q, + ^±[(^B-E)a] . (47) 
a dr 



The relation between and the general gravitational potential (j) obeys the Poisson 

equation [63| [6i]: 

= -k^cp = iTrGpa^6^^^. (48) 

where k is the (comoving) wavenumber of Fourier mode. As discussed above, the Poisson 
equation is valid only for scales small compared to the Hubble radius 1/H while on scales larger 
than the Hubble scale the growth of matter density perturbations is frozen. Hence, S^^^ can not 
be regarded as the observable matter density perturbation on scales comparable to the Hubble 
scale. Therefore, the observable density perturbation on both the small-scale and the large-scale 
modes can not be described directly by 6^^^ even though it is a gauge-invariant quantity. 

Contrary to 6^/^ , the other gauge invariant perturbation S^^^ has some important attractive 
features with respect to observability. These are summarized as follows: 

• It reduces to the Newtonian gauge perturbation 5^'-' ie it corresponds to a frame which 
respects the isotropic expansion of the universe and is therefore more appropriate for 
description of large scale perturbations. This reduction also simplifies the calculation of 
this perturbation. 

• It drives a scale dependent modification of the Poisson equation for the gauge invariant 
potential <j). Indeed, the time-time part of the linearized Einstein equation gives [65t [UP] 

y2 ^ _ ^a^a^ ^^^^ _^2^ _ ^a^a^ ^ ^ AnGpa^dg''' . (49) 
a a a a 

Thus, the anticipated scale dependence on scales comparable to the Hubble scale is picked 
up by the perturbation 6^^- 

• It is gauge invariant as anticipated for any observable quantity. 

These features make the gauge invariant Sg^^ and the Newtonian gauge to which it 
reduces, an attractive choice for making theoretical calculations to obtain the gravitational 
potential and the matter density perturbation that can be directly compared with observations 
on large scales. However, these theoretically obtained quantities need to also be corrected 
for bias, redshift distortions (due to gravitational potential and peculiar velocities), lensing 
magnification and volume distortion |62]. 

4. Conclusion 

I have reviewed the consistency of recent geometric cosmological data with the ACDM 
cosmological model. The standard candle SnIa datasets considered in this study are consistent 
with ACDM and with standard rulers at a level of 95.4% {2a) or less for certain prior values of 
the matter density r^om in the range ilom £ [0.25, 0.35]. It is interesting that the ranking sequence 
based on consistency with ACDM is practically identical with the corresponding ranking based 
on consistency with standard rulers even though the two criteria are completely independent. 
Thus, despite the improvement of standard ruler and standard candle data quality during 
the last decade the consistency of ACDM with data has not decreased despite the fact that 
ACDM is a simple, specific and well defined model which appears as a measure-zero point in all 
generalized models. On the contrary its consistency seems to be improving with time as new 
and more accurate data appear. For example, the Constitution SnIa dataset which is a very 
recent compilation with a drastic improvement on the crucial nearby SnIa sample, is also one of 
the most consistent datasets with both ACDM and standard rulers. 

Despite of its excellent consistency with both SnIa standard candles and CMB-BAO standard 
rulers, ACDM has to face potential challenges from other cosmological data [18j (e.g. large scale 
velocity flows, galaxy and cluster halo profiles, peculiar features of CMB maps etc.) which 
may lead the quest for the properties of dark energy to interesting surprises in the near future. 



Such surprises may also come from future standard candle observations or standard ruler CMB 
experiments (e.g. Planck [66]) which are expected to significantly improve the accuracy of the 
constraints discussed in the present study. 

I have also discussed dynamical probes of the cosmological constant and reviewed a scale 
dependent parametrization of the growth rate / = '^^^^^ that is free from the sub-Hubble 
approximation of the standard parametrization ()16p . This parametrization described by 
equation (|37p approximates very well the full linear general relativistic solution up to horizon 
scales and differs significantly from the standard parametrization on scales larger than about 
100h~^ Mpc. This parametrization describes well the growth of matter density perturbations on 
large scales in the Newtonian gauge. 

Even though the matter density perturbation is a gauge dependent quantity, the Newtonian 
gauge has certain advantages that make it particularly useful when making a direct comparison 
of theoretical predictions with observations. These advantages are summarized as follows: 

• The evolution of perturbations in the Newtonian gauge picks up the anticipated scale 
dependence when the Hubble scale is approached. 

• The time slicing of the Newtonian gauge respects the isotropic expansion of the background 
and therefore it is more appropriate to describe large scale fluctuations. On the other hand, 
the synchronous gauge corresponds to free falling observers at all points and therefore 
it is physically relevant on small scales where the two gauges give identical evolution of 
cosmological perturbations. 

• Gauge invariant generalizations of the density fluctuation 6m and the gravitational 
potential reduce to the corresponding Newtonian quantities in the Newtonian gauge. This 
identiflcation with gauge invariant quantities is a prerequisite for the direct observability of 
the density perturbation and the gravitational potential in the Newtonian gauge 

Even though additional corrections need to be made to the theoretically predicted value of the 
density perturbation evolution along the lines of Ref. [62] , it is clear that the Newtonian gauge 
offers a good starting point for the direct comparison of theoretical predictions of dark energy 
models with observations. 
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